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1.  Iterative  Substructuring  with  Lagrange  Multipliers  for 
Coupled  Fluid-Solid  Scattering 

Jan  Mandel1 

1.  Introduction.  In  [9],  we  have  proposed  an  iterative  method  for  the  solution  of 
linear  systems  arizing  from  finite  element  discretization  of  the  time  harmonic  acoustics 
of  coupled  fluid-solid  systems  in  fluid  pressure  and  solid  displacement  formulation. 
The  method  extended  the  FETI-H  method  for  the  Helmholtz  equation  [4,  6,  7,  12] 
to  coupled  fluid-elastic  acoustics.  In  this  paper,  we  investigate  a  stabilization  of  the 
discrete  coupled  system  for  the  case  when  the  solid  scatterer  is  at  resonance  and 
investigate  computationally  the  convergence  of  the  iterative  substructuring  method 
for  the  modified  system. 

The  main  idea  of  the  method  of  [9]  is  as  follows.  The  fluid  and  the  solid  domains  are 
decomposed  into  non-overlapping  subdomains.  Continuity  of  the  solution  is  enforced 
by  Lagrange  multipliers.  To  prevent  singular  or  nearly  singular  subdomain  matrices 
due  to  resonance,  the  continuity  conditions  betweeen  the  subdomains  are  replaced 
by  artificial  radiation-like  conditions.  Because  original  degrees  of  freedom  are  coupled 
across  the  wet  interface,  the  system  is  augmented  by  duplicating  the  degrees  of  freedom 
on  the  wet  interface  and  adding  equations  enforcing  the  equality  of  the  original  and  the 
duplicate  degrees  of  freedoom.  The  original  degrees  of  freedom  can  then  be  eliminated 
subdomain  by  subdomain  and  the  resulting  system  is  solved  by  Krylov  iterations 
preconditioned  by  a  Galerkin  correction  on  a  subspace  consisting  of  plane  waves  in  each 
subdomain.  In  each  iteration,  the  method  requires  the  solution  of  one  independent 
acoustic  problem  per  subdomain,  and  the  solution  of  a  coarse  problem  with  several 
degrees  of  freedom  per  subdomain.  The  number  of  iterations  in  was  most  cases  about 
the  same  as  the  number  of  iterations  of  the  FETI-H  method  for  the  related  Helmholtz 
problem  with  Neumann  boundary  condition  instead  of  an  elastic  scatterer,  which  was 
explained  by  numerical  decoupling  of  the  fluid  and  the  elastic  fields  in  the  stiff  scatterer 
limit. 

In  this  article,  we  propose  a  new  artificial  radiation-like  condition  on  the  wet 
interface,  and  we  observe  in  computational  tests  that  that  it  it  prevents  deterioration 
of  convergence  in  the  case  of  one  solid  subdomain  at  resonance.  We  also  investigate 
the  sensitivity  of  the  method  to  variants  of  artificial  radiation  condition  between  the 
elastic  subdomains. 

Our  radiation- like  condition  between  elastic  subdomains  has  been  inspired  by  [2], 
which  generalized  the  alternating  method  of  [5]  to  elasticity.  Iterative  methods  con¬ 
sisting  of  alternating  solution  in  the  fluid  and  the  solid  region  are  known  [1,  3].  In  [3], 
the  alternating  method  of  [5]  was  extended  to  the  coupled  problem,  with  the  wet  inter¬ 
face  conditions  replaced  by  their  complex  linear  combinations.  The  resulting  iterative 
algorithm  needs  either  access  to  normal  derivatives  or  additional  variables  on  the  wet 
interface.  Our  radiation-like  condition  on  the  wet  interface  is  obtained  by  a  simple 
modification  of  the  coupled  system  matrix,  resulting  in  an  equivalent  algebraic  system. 
Since  this  process  is  unrelated  to  the  substructuring  method  at  hand,  it  may  be  of 
independent  interest. 
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2.  The  scattering  problem.  We  need  to  describe  the  scattering  problem  and 
the  discretization  used.  This  material  is  standard  [11,  13]  and  it  is  included  only  for 
completeness  and  to  introduce  the  notation. 

We  consider  an  acoustic  scattering  problem  with  an  elastic  scatterer  completely 
immersed  in  a  fluid.  Let  f 1  and  fie  be  bounded  domains  in  SR",  =  2,3,  f 1e  C  f 2,  and 
let  0/  =  f2  \  fie,  cf.,  Figure  5.1.  Let  v  denote  the  exterior  normal  of  fie.  Let  3fi  be 
decomposed  into  disjoint  subsets,  3fi  =  U  T„  U  ro.  The  domain  fi j  is  filled  with  a 
fluid.  The  acoustic  pressure  at  time  t  is  assumed  to  be  of  the  form  R epelwt,  where  p  is 
complex  amplitude  independent  of  t.  The  amplitude  p  is  governed  by  the  Helmholtz 
equation 

Ap  +  k2p  =  0  in  fi  f,  (2.1) 

with  the  boundary  conditions 

dv  dv 

P  =  Po  on  Tj,  —  =0onr„,  +  ikp  =  0  on  T0,  (2.2) 

where  k  =  cj/c/  is  the  wave  number  and  c/  is  the  speed  of  sound  in  the  fluid. 
The  boundary  conditions  (2.2)  model  excitation,  sound  hard  boundary,  and  outgo¬ 
ing  boundary,  respectively.  The  amplitude  of  the  displacement  u  of  the  elastic  body 
occupying  the  domain  fie  satisfies  the  elastodynamic  equation 

V  •  r  4-  u>2peu  =  0  in  fie,  (2.3) 

where  r  is  the  stress  tensor  and  pe  is  the  density  of  the  solid.  For  simplicity,  we 
consider  an  isotropic  homogeneous  material  with 

r  =  \I{V -u)+2pe{u),  e^u)  =  (2-4) 

where  A  and  p  are  the  Lame  coefficients  of  the  solid. 

Let  T  =  9fie  be  the  wet  interface.  On  T,  the  fluid  pressure  and  the  solid  displace¬ 
ment  satisfy 

V  ■  U  =  — — V  ■  T  ■  V  =  —p,  VXT'V=  0,  (2.5) 

Pfioz  av 

where  pf  is  the  fluid  density.  The  interface  conditions  (2.5)  model  the  continuity 
of  normal  displacement,  the  balance  of  normal  forces,  and  zero  tangential  tension, 
respectively. 

We  use  the  following  variational  form.  Define  the  spaces  Vf  =  {q  €  iL1(fi/)  |  q  = 
0  on  Tj},  Ve  =  (iL 1  (fie))",  where  Hl  is  the  Sobolev  space  of  generalized  functions 
with  square  integrable  generalized  first  derivatives.  Assuming  that  po  on  Tj  is  ex¬ 
tended  to  a  function  in  U1(fi/),  multiplying  equation  (2.1)  by  a  test  function  q  £  Vf, 
equation  (2.3)  by  a  test  function  u  £  Ve,  and  integrating  by  parts,  we  obtain  the 
following  variational  form  of  (2.1)  -  (2.5):  Find  p  such  that  p  —  po  £  Vf,  and  u  £  Ve 
such  that  for  all  q  £  Vf  and  all  v  £  Ve, 


■  J  VpVq  +  k2  jpq  —  ik  /  pq-u 2  /"'<  v-u)q 


r  a 


J  (A (V  •  u)(V  •  v)  +  2 pe{u)  :  e(v))  +  J2  J  peu  ■  v.  -  J p(v  ■  v) 


r 


=  0, 

=  0. 
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We  replace  V;  and  Ve  with  conforming  finite  element  spaces  and  obtain  the  algebraic 
system 


— K  j  +  fc2M  f  —  ikGf 

—pfOJ2  T 

P 

r 

— T' 

— Ke  +  c j^lVIg 

u 

0 

In  the  coupled  system  (2.6),  p  and  u  are  the  vectors  of  the  (values  of)  degrees 
of  freedom  of  p  and  u,  i.e.,  p  and  u  are  the  finite  element  interpolations  of  p  and  u, 
respectively.  The  matrix  blocks  in  (2.6)  are  defined  by 


p'K/q  =  J  Np-Nq,  p'M/q  =  J  pq, 

£lf  £2  f 

p  G/q  =  J pq,  u'Kev  =  J  (A(V  ■  w)(V  ■  v)  +  2pe(u)  :  e(v)), 
ra  ne 

u'Mev  =  J  pe(u-v),  p'Tv  =  Jp(v-v). 


3.  Iterative  Substructuring.  In  this  section,  we  summarize  the  iterative  method 
following  [9].  Further  details  and  a  development  of  the  method  starting  from  FETI-H 
can  be  found  in  [9].  The  present  method  differs  in  the  more  general  choice  of  artificial 
radiation  condition  between  elastic  subdomains. 

The  fluid  and  solid  domains  are  decomposed  into  nonoverlapping  subdomains  that 
consist  of  unions  of  elements, 


Nf  _  _  Ne  _ 

%  =  1J  K  ne  =  (j  K  (3.i) 

s—l  s=  1 

The  fields  and  vectors  of  degrees  of  freedom  corresponding  to  and  fl*  are  denoted 
by  ps,  us,  ps  and  us,  respectively.  The  normal  vector  to  <9fis  is  denoted  by  vs . 

The  Helmholtz  equation  (2.1)  is  then  equivalent  to  the  same  equation  in  each  of 
the  subdomains  Uj,  with  the  interface  conditions 

pS=pt’  +  ^  =  on50/n50/-  O-2) 

Similarly,  the  elastodynamic  equation  (2.3)  is  equivalent  to  the  same  equation  in  each 
of  the  subdomains  0* ,  with  the  continuity  of  the  displacement  and  the  traction  on  the 
intersubdomain  interfaces, 

us=ut,  t(us)vs  4-  T(ut)ut  =  0,  on  ndQ*.  (3.3) 

The  continuity  of  the  pressure  and  the  displacement  will  be  enforced  by  Lagrange 
multipliers. 
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Define  subdomain  matrices  by  subassembly, 


us,Mfv 


=  f  Vp- 

J 

Vg, 

Qsf 

=  / 

PQ, 

9f2“nra 

II 

■ v ) 

5,M^qs  =  J  pq, 


us,K*v: 


=  J  A (V  •  u)(y  ■  v)  +  2 fie(u)  :  e(v), 


pr,Trsvs  =  J  p(v  v). 
We  will  use  vectors  consisting  of  all  subdomain  degrees  of  freedom, 


‘  P1  ‘ 

■  u1  ■ 

15 

,  U  = 

C 

£  " 

w 

I—1 

0  ■ 

'  K) 

0 

,  Ke  =  diag(K®)  = 

o  . 

W 

0  • 

“5? 

I _ 

and  the  corresponding  partitioned  matrices, 

K  /  =  diag(Kj)  = 

The  matrices  M/,  G /,  and  Me  are  defined  similarly,  and 
T  =  (Trs)r.s  =  ' 

I 

rj^Nf,l  rj^Nf,Ne 


Let  N /  and  Ne  be  the  matrices  with  0,1  entries  of  the  global  to  local  maps  corre¬ 
sponding  to  the  decompositions  of  £lf  and  fle,  respectively,  cf.,  (3.1),  so  that 

K/=N^K/N/,  Ke  =  N'eKeNe. 

Let  B/  =  (Bj, . . . ,  B^)  and  Be  =  (Bj, . . . ,  B^e)  be  matrices  of  full  rank  such  that 
the  conditions  B^p  =  0  and  Beu  =  0  express  the  constraint  that  the  values  of  the 
same  degrees  of  freedom  on  two  different  subdomains  coincide,  that  is,  B^p  =  0  ^=>- 
p  =  N/p  for  some  p,  and  Beu  =  0  •<=>•  u  =  Neu  for  some  u.  See  [8]  for  details 
on  the  construction  of  such  matrices  with  entries  0,  ±1.  Here,  we  use  the  matrices 
from  [8]  and  orthogonalize  their  rows  for  numerical  stability;  the  resulting  matrices 
are  still  sparse. 

Multiplying  the  second  equation  in  (2.6)  by  w2/)/  to  symmetrize  the  off-diagonal 
block  and  introducing  Lagrange  multipliers  A /  and  Ae  for  the  constraints  B/p  =  0 
and  Beu  =  0,  we  get  the  system  of  linear  equations  in  block  form, 


'  -Kf  +  k2Mf  -  ikG 

-cj2pfT 

ffl 

0 

P 

f 

-w2PiT' 

co2Pf(-Ke  +CJ2  Me) 

0 

B'e 

u 

0 

B  / 

0 

0 

0 

A/ 

0 

0 

Be 

0 

0 

.  . 

0 

(3.4) 
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where  N'r  =  r.  Similarly  as  in  [8],  it  can  be  shown  that  the  system  (3.4)  is  equivalent 
to  (2.6)  in  the  sense  that  (p,  u)  is  a  solution  of  (2.6)  if  and  only  if  (p,  u,  A/,  Ae)  with 
p  =  N/p,  u  =  Neu,  is  a  solution  of  (3.4)  for  some  A /  and  Ae. 

Using  the  properties  of  the  global  to  local  maps  N/  and  Ne,  it  is  easy  to  see  that 
(p,  u,  A/,  Ae)  is  a  solution  of  (3.4)  if  and  only  if  p  =  N/p  and  u  =  Neu,  where  (p,  u) 
solves  (2.6). 

We  will  want  to  eventually  eliminate  the  variables  p  and  u.  But  the  matrices 
— K f  +  k2M.f  and  — Ke  +  w2 Me  may  be  singular  due  to  resonance.  For  this  reason, 
the  continuity  of  normal  derivative  and  traction  between  subdomains  are  replaced  by 
artificial  radiation  conditions, 

ps  +  iastk-^—  =ut+  iatsk-^-r  on  <9fK  n  (3.5) 

dvs  duf  ;  ;  v  ' 

and 

us  +  iatsaT(us)vs  =  u*  +  iast ar{ui:)vi: ,  on  n  90*-  (3-6) 

Here,  ast  =  ±1  or  0,  ast  =  —  ats ,  and 

a  =  a0u>\/pe(\  +  2  n).  (3.7) 

If  ast  =  ±1,  the  interface  condition  (3.5)  allows  a  plane  wave  to  pass  in  one  nor¬ 
mal  direction  through  the  interface  between  the  subdomains.  Similarly,  by  a  simple 
computation,  the  condition  (3.6)  with  «o  =  1  is  satisfied  by  the  plane  pressure  wave 


u(x)=decpd  X,  |d|  =  1,  cp  =  \  (3-8) 

V  Pe 

in  one  of  the  normal  directions,  d  =  ±vs.  An  alternative  form  of  (3.6)  is 

vs (us  ■  vs)  +  i<jts oit(us)vs  =  vt(ut  ■  vt)  +  iasi lctT{ut)vt  on  cfflj  n  dillj ,  (3.9) 

which,  for  a  from  (3.7)  with  «o  =  1,  is  also  satisfied  by  the  pressure  wave  (3.8)  in 
normal  direction.  In  [9],  the  condition  (3.9)  with  a  =  upe  was  used. 

This  change  of  intersubdomain  interface  conditions  corresponds  to  replacing  the 
subdomain  matrices  — K f  +  k2M.f  and  — Ke  +w2Me  by  regularized  matrices 

Aj  =  —  Kj  +  k2NIf  +  ikGf  +  R/, 

Ae  =  — Ke  +  w2Me  +  Re, 


where  the  regularization  matrices  are  given  by 


R/  =  diag(Rj), 


between  fluid  subdomains,  and 


Re  =  diag(Rg), 


ive  , 

us,RgVs  =  ia  °st  / 

■4- - 1  J 


u  ■  V , 


t= 1 
t^s 


dQf  (IdQ* 


(3.10) 
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between  elastic  subdomains  for  the  interface  condition  (3.6)  and  by 

Ne  . 

Re  =  diag(R*),  us,R*vs  =  ia  ^  ast  /  {vs  ■  u){vs  ■  v),  (3.11) 

51 

if  (3.9)  is  used. 

It  is  shown  in  [6]  for  the  Helmholtz  equation  that  if  for  a  given  s,  all  ast  >  0  or 
all  ust  <  0  with  some  ast  0,  then  is  invertible.  The  case  of  elastic  subdomains 
is  similar.  For  details  on  strategies  for  choosing  ast  to  guarantee  this,  see  [6].  In  our 
computations,  we  simply  choose  ast  =  +1  if  s  >  t,  as^  =  —  1  if  s  <  t. 

Because 

N| Ry Ny  =  0,  N' ReNe  =  0, 

the  effect  of  adding  the  matrices  R /,  R/  cancels  in  the  assembled  system,  and  the 
system  (3.4)  is  equivalent  to 

Af  —LJ2pfT  Bj  0  p  r 

-cj2pff'  u2pfAe  0  B'e  u  _  0  .  . 

B  f  0  0  0  A  f  0 

0  Be  0  0  J  L  Ae  J  L  0  . 

Eliminating  the  original  degrees  of  freedom  at  this  point  does  not  result  in  inde¬ 
pendent  computation  in  each  subdomain,  because  of  coupling  of  degrees  of  freedom 
across  the  wet  interface  by  the  matrix  T.  Hence,  we  first  duplicate  the  interface  de¬ 
grees  of  freedom  as  follows.  Since  the  value  of  Tu  depends  on  the  values  of  u  on  T 
only,  we  have 

Tu  =  TJeur,  ur  =  JeU, 

where  Je  is  the  matrix  of  the  operator  of  embedding  a  subvector  that  corresponds  to 
degrees  of  freedom  on  T  into  u  by  adding  zero  entries.  Similarly, 

T'p  =  T'J/pr,  Pr  =  J/P- 

Therefore,  we  obtain  the  augmented  system  equivalent  to  (3.12), 


A/ 

0 

B  / 

0 

0 

~UI2pf  TJe 

0 

u2pfAe 

0 

B'e 

— ui2  pfT'J  f 

0 

B  / 

0 

0 

0 

0 

0 

0 

Be 

0 

0 

0 

0 

J/ 

0 

0 

0 

-I 

0 

0 

J'e 

0 

0 

0 

-I 

Because  the  variables  in  a  coupled  system  typically  have  vastly  different  scales,  we 
use  symmetric  diagonal  scaling  to  get  the  scaled  system 


A/ 

0 

B '/ 

0 

o 

o 

Ae 

0 

B' 

— T'J  / 

B  / 

0 

0 

0 

0 

o 

Be 

0 

0 

0 

J  / 

0 

0 

0 

-I 

0 

J'e 

0 

0 

0 

TJe  p  f 

0  u  0 

0  A/  =  0 

0  Ae  O’ 

0  Pr  0 

-I  J  L  ur  J  L  0  . 


(3.14) 
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where  the  matrices  and  the  vectors  scale  as  A /  =  D/A/D/,  Ae  =  cj2pf DeAeDe, 
T  =  u2pfDfTDe,  B/  =  E/B/D/,  Be  =  EeBeDe,  f  =  D/r,  p  =  D/p,  u  =  Deu, 
Xf  =  D/A/,  Ae  =  DeAe.  The  scaling  matrices  D/,  De,  E /,  and  Ee,  are  diagonal. 
We  have  chosen  scaling  matrices  with  positive  diagonal  entries  such  that  the  absolute 
values  of  the  diagonal  entries  of  A/  and  Ae  are  one  and  the  I2  norms  of  the  columns 
of  Be  and  B/  are  one. 

Computing  p  and  u  from  the  first  two  equations  in  (3.14)  gives 

p  =  Ay  x(r  —  B/A/  +  TJeur)  (3.15) 

u  =  A“1(— B(,Ae  +T,J/pr)  (3.16) 

Substituting  p  and  u  from  (3.15),  (3.16)  into  the  rest  of  the  equations  in  (3.14),  we 
obtain  the  reduced  system 

Fx  =  b,  (3.17) 

where 

B/A^B).  0  0  -B/A^TJe 

0  Be  A~1Be  -BeAj1T'J  / 

—  j'yAy^B'y  0  ~I  J'yA^TJe 

0  -3ek~lBe  JeA-iT'J  /  -I 

and 

B/A^f 
0 

0 

In  equation  (3.18),  the  first  diagonal  block  B/A^B'y  is  exactly  same  as  in  the 
FETI-H  method  for  the  Helmholtz  equation.  The  second  diagonal  block  B/Ay  1By  is 
the  analogue  of  FETI-H  for  the  elastodynamic  problem. 

Evaluating  the  matrix  vector  product  Fx  requires  the  solution  of  one  independent 
problem  per  subdomain,  because 

’  A/  1  [  -B/q  1 

F  Xe  =  ,  where  H  =  A/11(“B/A/+TJefir)’ 

Pr  J'/q-Pr  l  v  =  A71(-B'eAe+T'J/pr). 

ur  J'ev  -  ur 

The  iterative  method  then  consists  of  solving  the  linear  system  (3.17)  by  GMRES 
preconditioned  by  a  subspace  correction  as  follows.  Let  Q  be  a  matrix  with  the  same 
number  of  rows  as  F  and  linearly  independent  columns.  The  columns  of  Q  form  the 
basis  of  the  coarse  space.  The  orthogonality  condition 

Q'(Fx  -  b)  =  0,  (3.19) 

is  enforced  through  the  iterations  by  adding  a  correction  from  the  coarse  space  in  each 
iteration.  That  is,  GMRES  is  applied  to  the  preconditioned  system 


PFx  =  Pb 


(3.20) 
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where  P  =  (I  -  Q(Q'FQ)_1Q'F)  and  the  initial  approximation  x  =  Q(Q'FQ)-1b 
satisfies  (3.19).  Because  the  increments  are  in  the  range  of  P  and  Q'FP  =  0,  all 
iterates  satisfy  (3.19). 

We  choose  the  matrix  Q  of  the  form 


'  D/B/diag(Y^)s 

0 

0 

0 

0 

DeBediag(Y|)s 

0 

0 

0 

0 

D/J'/diag(Zps 

0 

0 

0 

0 

DeJ'ediag(Z®) 

(3.21) 


and  orthogonalize  its  columns  by  the  QR  algorithm.  For  a  fluid  subdomain  ilj. 
we  choose  as  the  matrix  of  columns  that  are  discrete  representations  of  plane 
waves  in  a  small  number  of  equally  distributed  directions,  or  discrete  representation 
of  the  constant  function.  For  a  solid  subdomain  fi®,  the  columns  of  Yf  are  discrete 
representations  of  plane  pressure  and  shear  waves,  or  of  rigid  body  motions.  The 
matrices  Z)-  and  Z®  are  chosen  in  the  same  way  as  Y)  and  Y®.  See  [9]  for  further 
details  and  a  discussion  why  the  method  for  the  coupled  problem  can  be  expected  to 
perform  about  as  FETI-H  for  the  fluid  and  the  elastic  parts  separately. 

4.  Radiation-Like  Condition  on  the  Wet  Interface.  For  some  frequencies, 
the  matrix  — Ke  +  lj2 Me  in  the  coupled  system  (2.6)  will  be  singular.  The  inverse  of 
this  matrix  is  required  by  the  method  if  there  is  only  one  elastic  subdomain.  Therefore, 
in  this  case,  we  replace  (2.6)  by  an  equivalent  system,  obtained  by  adding  to  the  second 
block  of  equations  a  linear  combination  of  the  first  block  in  such  a  way  that  the  term 
added  to  — Ke  +  cu2Me  resembles  a  radiation  condition: 


— Kj  +  fc2Mj  —  ikGf  —pfLU2  T 

-T'  +*r^T'(-K/  +k2Mf  -ikGf)  -Ke  +  w2Me  +  i/3T'T 

where 

r  r  1 


To  obtain  an  added  term  with  consistent  physical  units  and  similar  to  the  artificial 
radiation  condition  (3.10),  we  choose 

(4.2) 

In  the  case  of  more  than  several  fluid  subdomains  and  one  elastic  subdomain,  this 
process  is  easily  implemented  using  the  local  subdomain  matrices  for  T.  For  more 
than  one  elastic  subdomain,  computational  experiments  indicate  that  introducing  an 
artificial  radiation  condition  on  the  wet  interface  is  not  necessary. 

5.  Computational  results.  Computational  results  showing  scalability  of  the 
method  were  presented  in  [9].  Here,  we  focus  on  the  performance  of  the  method  when 
k  equals  or  is  close  to  a  value  that  makes  some  of  the  subdomain  matrices  singular, 
and  for  different  choices  of  the  radiation-like  condition  between  elastic  subdomains. 
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Figure  5.1:  Model  2D  Problem 

r„ 


We  consider  a  model  2D  problem  with  a  scatterer  in  the  center  a  waveguide,  cf., 
Fig.  5.1.  The  fluid  domain  Llf  is  a  square  with  side  1  m,  filled  with  water  with  density 
Pf  =  1000  kg  m~3  and  speed  of  sound  Cf  =  1500ms-1.  The  scatterer  is  a  square  in 
the  center  of  the  fluid  domain,  consisting  of  aluminum  with  density  pe  =  2700  kgm~3 
and  Lame  elasticity  coefficients  A  =  5.5263.10 10  Nm~2,  p  =  2.595.10 10  Nm~2.  The 
domain  is  discretized  with  a  mesh  of  200  by  200  bilinear  elements.  The  coarse  space 
consists  of  8  plane  waves  in  the  fluid  subdomains  and  4  plane  pressure  waves  and 
4  plane  shear  waves  in  the  solid  subdomains  (the  first  two  blocks  of  the  matrix  Q 
in  (3.21)).  The  same  number  of  coarse  space  functions  is  used  for  the  coarse  space  for 
the  wet  interface  (the  last  two  blocks  of  the  matrix  Q).  The  iterations  are  terminated 
when  the  relative  residual  reached  10-6.  Then  the  scaled  residual  in  the  original 
variables  are  checked, 


|dj  KjjZj 

Res  =  max  — r 

1  2-^j  \zj\ 


(5.1) 


where  z,  K,  and  d  are  the  solution  vector,  the  matrix,  and  the  right-hand-side,  re¬ 
spectively,  of  the  coupled  system  (2.6).  In  all  cases  when  the  iterations  converged, 
this  scaled  residual  was  of  the  order  10-6  to  10-7. 

Figure  5.2  shows  the  number  of  iterations  for  varying  constant  «o  in  the  artificial 
radiation  condition  (3.6)  between  elastic  subdomains.  The  scatter  was  size  0.4  by  0.4 
and  the  fluid  and  the  elastic  domains  were  decomposed  into  4  subdomains  each.  One 
can  see  that  the  number  of  iterations  for  a0  =  1  is  slightly  larger  over  all  frequen¬ 
cies,  while  the  iterations  diverge  for  frequencies  equal  to  or  very  close  the  resonance 
frequencies  for  «o  <  10-4. 

The  number  of  iterations  for  the  same  test  problem  and  the  artifical  radiation 
condition  (3.9)  was  almost  exactly  same  (not  shown). 

Figure  5.3  reports  the  number  of  iterations  for  the  same  problem  with  the  artificial 
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Figure  5.4:  Number  of  iterations  for  different  80  in  artificial  radiation  on  wet  interface 
by  (4.1)  and  (4.2) 


radiation  condition  (3.6)  between  elastic  subdomains,  but  instead  of  orthogonalization 
of  the  rows  of  the  matrices  B  and  the  columns  of  the  matrix  Q,  bases  are  selected 
as  linearly  independent  subsets.  There  are  more  iterations  required  and  divergence 
occurs  for  more  frequencies  and  larger  values  of  the  parameter  «o- 

Figure  5.4  reports  the  number  of  iterations  for  decreasing  strength  80  of  the  artifi¬ 
cial  radiation-like  term  on  the  wet  interface.  The  scatterer  was  size  0.2  by  0.2,  forming 
one  elastic  subdomain,  and  the  fluid  domain  was  decomposed  along  the  midlines  of  the 
square  into  4  subdomains.  One  can  see  that  the  choice  80  =  1  increases  the  number  of 
iterations  significantly  over  all  frequencies,  while  for  80  —  10-6,  the  iterations  diverge 
for  frequencies  equal  to  or  close  to  a  resonance  frequency.  The  elastic  subdomain 
in  this  experiment  is  of  the  same  size  as  the  elastic  subdomain  for  the  examples  in 
Figure  5.2,  so  the  resonance  frequencies  are  same. 
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